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ABSTRACT 

We present images of the jets in the nearby radio galaxy NGC315 made with the 
VLA at five frequencies between 1.365 and 5 GHz with resolutions between 1.5 and 
45arcsec FWHM. Within 15arcsec of the nucleus, the spectral index of the jets is 
a — 0.61. Further from the nucleus, the spectrum is flatter, with significant transverse 
structure. Between 15 and 70 arcsec from the nucleus, the spectral index varies from 
«0.55 on-axis to «0.44 at the edge. This spectral structure suggests a change of 
dominant particle acceleration mechanism with distance from the nucleus and the 
transverse gradient may be associated with shear in the jet velocity field. Further from 
the nucleus, the spectral index has a constant value of 0.47. We derive the distribution 
of Faraday rotation over the inner ±400 arcsec of the radio source and show that 
it has three components: a constant term, a linear gradient (both probably due to 
our Galaxy) and residual fluctuations at the level of 1 - 2 radm" 2 . These residual 
fluctuations are smaller in the brighter (approaching) jet, consistent with the idea that 
they are produced by magnetic fields in a halo of hot plasma that surrounds the radio 
source. We model this halo, deriving a core radius of «225 arcsec and constraining its 
central density and magnetic-field strength. We also image the apparent magnetic-field 
structure over the first ±200 arcsec from the nucleus. 

Key words: galaxies: jets - radio continuum:galaxies - galaxies: individual: NGC 315 
- magnetic fields - polarization - MHD 



1 INTRODUCTION 



The giant FRI (iFanaroff fc Rilevl Il974) r adio source 
NGC 315 was first imaged bv iBridle et al.1 <1976ft . who 

showed that it has an angular size of nearly 1°. The 
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IVenturi et^ll993t ICotton et al]ll999l: IXu et alJl200* . X- 

ray emissio n from the first 10 arcsec of the main je t was 
detected by IWorrall. Birkinshaw fc Hardcastlel i2003^) , but 
no optical emission from this region has yet been reported. 

The source is associ ated with a giant elliptical galaxy 
at a redshift of 0.01648 llTraeer et alJ l2000). giving a scale 
of 0.335 kpc/arcsec for our adopted cosmology (Hubble 
constant Ho = 70kms _1 Mpc -1 , Qa = 0.7 and Qm = 



0.3). NGC 3 15 is a member of a group or poo r cluster 
of galaxies dNoltheniusI Il993t iMiller et, al.1 1200* located 
in o ne of the filaments of the Pisces-Pers eus superclus- 
ter jEnfilin et alJ 1200 it iHuchra et al.1 Il99* . HST images 
ijVerdoes^^rmet alJll99* show a 2.5-arcsec diameter dust 
lane and a nuclear point source. The dust lane is associ- 
ated with a disk of ionized gas w hich is probably in ordered 
rotation ijNoel-Storr et alJl200* . CO emission, also with a 
line p rofile indicating rotation, was detected bv lLeon et alJ 
(2003). The inferred mass of molecular hydrogen is (3.0 ± 
0.3) x 10 8 M Q and the cold gas is likely to be cospatial with 
th e dust. HI absorption ag ainst the nucleus was detected 
bv lvan Gorkom et alJ il989T) . There is evidence fo r a weak 
polar i zed broad Ha line in the nuc l ear spectrum dHo et al 



I l997t iBarth. Filippenko fc Moranl Il999t iNoel-Storr et al 
2003). Hot gas associated wit h the galaxy has been imaged 
using ROSAT and Chandra, JWqrrall fc Birkinshawl lioool 
IWorrall. Birkinshaw fc Hardcastldl2003ft . 



E-mail: rlaing@eso.org 



Within ~ 90 arcsec of the nucleus, the jets in NGC 315 
are initially narrow, then expand rapidly ("flare") and re- 



2 R.A. Laing et al. 



collimate feridlelll98i ICanvin et alJl2005l) . We have mod- 
elled the inner ±70arcsec of this flaring region as a two- 
sided, symmetrical, relativistic flow, fitting to deep, high- 
resolution VLA observations at 5 GHz in order to derive the 
three-dimensional distributions o f velocity, proper em issiv- 
ity and magnetic-field structure (ICanvin et al .1 12005). Our 
main conclusions are as follows. 

(i) The jets are inclined by 38° ± 2° to the line of sight. 

(ii) Where they first brighten, their on-axis velocity is 
(5 = v/c ~ 0.9. They decelerate to /3 ~ 0.4 between 8 and 
18kpc from the nucleus (15 - 33arcsec in projection) and 
the velocity thereafter remains constant. 

(iii) The ratio of the speed at the edge of the jet to its 
value on-axis ranges from ~0.8 close to the nucleus to ~0.6 
further out. 

(iv) The longitudinal profile of proper emissivity is split 
into three power-law regions separated by shorter transition 
zones and the emission is intrinsically centre-brightened. 

(v) To a first approximation, the magnetic field evolves 
from a mixture of longitudinal and toroidal components to 
predominantly toroidal by 26kpc (48arcsec in projection). 

(vi) Simple adiabatic models fail to fit the emissivity vari- 
ations. 

In the present paper, we investigate the energy spec- 
trum of the relativistic particles in th e jets of NGC 315 i n 
the context of the models developed bv lCanvin et al.l (2005). 
We use VLA observations at frequencies between 1.365 and 
5 GHz 1 to derive the spectrum of the jets at resolutions of 
5.5 and 1.5 arcsec and relate the observed spectral gradients 
to velocity, emissivity and field structure. A separate paper 
(Worrall et al., in preparation) will describe the radio struc- 
ture of the main jet at high resolution and its relation to 
new Chandra images. 

We also determine the variations of Faraday rotation 
over the jets and test the hypothesis that these result from 
magnetic-field irregularities in hot, X-ray emitting plasma 
associated with the surrounding group of galaxies. Finally, 
we determine the apparent magnetic-field str ucture of the 
jets o n scales larger than those covered by ICanvin et al.l 
l|2005h . 

In Section |5J we describe the observations and their 
reduction. The total-intensity images are presented in Sec- 
tion 03 and we use them to derive distributions of spectral 
index in Section[f] We then discuss the distributions of Fara- 
day rotation (Section and apparent magnetic-field struc- 
ture (Section derived from observations of linear polar- 
ization. Section |7| summarizes our main results. 



2 OBSERVATIONS AND IMAGES 
2.1 Observations 

VLA data were obtained at 4.985 GHz in the B , C/D and 
A/D configurations as described bv lVenturi et al.l l)l993f) and 
ICotton et ail Jl999h . These were supplemented by additional 
observations in the A and C configurations with a centre fre- 
quency of 4.860 GHz to give complete coverage of the spatial 

1 The 5-GHz observations are those discussed by ICanvin et alJ 
i2005l) 



Table 1. Record of VLA observations, v and Ai/ are the cen- 
tre frequency and bandwidth, respectively, and t is the on-source 
integration time. 
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2001 
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scales accessible to the VLA in a single pointing. In order 
to map Faraday rotation, we observed at 1.365, 1.413, 1.485 
and 1.665 GHz in the B and C configurations of the VLA 
using a lower bandwidth. We also extracted observations 
in A configuration at 1.413 GHz from the VLA archive. 2 A 
journal of observations is given in Tabled 



2.2 Data reduction 

All of the data reduction was done in the MPS package. 
Initial amplitude and phase calibration were applied using 
standard methods and the flux-density scales were set us- 
ing observations of 3C 48 or 3C 286. Standard instrumen- 
tal polarization calibration was applied and the zero-points 
of E-vector position angle were determined using observa- 
tions of 3C 138 or 3C 286. The data for each configuration 
were first adjusted to a common phase centre in J2000 co- 
ordinates, imaged and self-calibrated separately. They were 
then concatenated in turn, starting with the widest config- 
uration. The slight difference in centre frequency between 
the datasets at 4.86 and 4.985 GHz was ignored (we show 
in Section [5.11 that this has a negligible effect on the anal- 
ysis of polarization) and we will refer to the combination 
as the "5 GHz dataset". At this frequency, the core showed 
signifi cant variability between observations (cf. lLazio et alJ 
2001) and the flux density of the unresolved component in 
the larger-configuration dataset was adjusted to match that 
observed with the smaller configuration when both were 
imaged at matched resolution. No core variability was de- 
tected at lower frequencies. A further iteration of phase 
self-calibration was done after each combination. Our final 
datasets are listed in Table £3 together with the minimum 
and maximum spatial scales they sample. 



2 In additio n, we re-analysed th e C-configuration dataset at 
8.4 GHz from Vent uri et al.Nl993l) . but poor weather during the 
observations precluded accurate absolute flux and polarization 
calibration, so we do not discuss them here. 
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Figure 1. Grey-scales of total intensity at low resolution, (a) 
WSRT image at 327 M Hz, showing the fu ll extent of the radio 
emission from NGC315 iMack et al.lll997T) . The restoring beam 
is 110 X 55arcsec 2 FWHM in position angle and the grey- 
scale range is - 75mjy (beam area) -1 . The box marks the area 
covered by panel (b) and Fig. |3 (b) The mean of all four VLA 
L-band images, tapered to a resolution of 45 arcsec FWHM and 
corrected for attenuation by the primary beam before averaging. 
The grey-scale range, — 5mjy (beam area)" 1 , has been chosen 
to emphasise the diffuse emission. In both panels, the cross marks 
the position of the core. 



Table 2. Final uv datasets. The columns are: (1) centre fre- 
quency, (2) array configurations used, (3) minimum and (4) max- 
imum spatial scales. 



Our observations were designed to give the maximum 
sensitivity for the inner jets of NGC 315 and were therefore 
taken with the pointing centre on or near the nucleus. As 
can be seen from Table H the maximum spatial scale sam- 
pled adequately at 5 GHz is ~300 arcsec. Even at L-band, 
the lobe associated with the counter-jet llMack et al]ll9971 
and Fig. is severely attenuated by the primary-beam re- 
sponse of the VLA and is not visible on our images. We did 
not recover the total flux density of the source at any fre- 
quency, so we estimated appropriate zero-spacing flux den- 
sities from the shortest-spacing visibility amplitudes. We 
made images at five resolutions: 45, 5.5, 2.35, 1.5 and 0.4 
arcsec FWHM, using similar baseline ranges at all frequen- 
cies and weighting the data in the uv plane as required. 
After imaging, we made both CLEAN and maximum-entropy 
deconvolutions. Although the latter algorithm gave slightly 
smoother images, it introduced a significant large-scale rip- 
ple parallel to the jet axis, whereas CLEAN gave a flat back- 
ground. We therefore show the CLEAN images, although we 
quote quantitative results only where the two deconvolu- 
tion methods agree. We also compared the I images made 
with and without zero-spacing flux densities and before and 
after subtraction of a local zero-level. None of these differ- 
ences led to significant changes in spectral index or degree of 
polarization compared with the errors quoted below. After 
deconvolution, all of the images were corrected for primary 
beam attenuation. We then took averages of the I images at 
1.365 - 1.665 MHz ("mean L-band images"). 

Data in Stokes Q and U were imaged without zero- 
spacing flux den sities and CLEANed. A first- order correction 
for Ricean bias dWardle fc Kronberdll97l was applied to 
the images of polarized intensity P = (Q 2 + U 2 ) 1 ^ 2 used to 
derive the degree of polarization p — P/I. 

The off-source noise levels at the centre of the field for 
the final images are given in Table (the 0.4-arcsec im- 
age at 5 GHz is discussed in detail by Worrall et al., in 
preparation, and is therefore not considered further here). 
Note that the wide-field L-band images at a resolution of 
5.5 arcsec are significantly affected by bandwidth smearing 
in their outer regions, images of point sources being broad- 
ened by a factor of 2 in the radi al direction at a distance of 
22 ar cmin from the phase centre dTavfor. Ulvestad fc Perfevl 
2004). This limitation needs to be taken into account only 
in the discussion of the source morphology on large scales 
(Section [3J. Measurements of spectra are restricted to the 
inner 200 arcsec of the field, where the effects of bandwidth 
smearing are <3% in peak intensity or image size for any 
of our frequency /resolution combinations. Our estimates of 
Faraday rotation, which extend to larger scales, should not 
be affected systematically by bandwidth smearing. 
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3 TOTAL INTENSITY 

A grey-sc ale of the large-sca le radio structure of NGC 315 at 
327 MHz jMack et al.lll997f) is given in Fig-EJa). Our obser- 
vations are sensitive only to emission from the region indi- 
cated by the box in this figure and our mean L-band image 
of this region at a resolution of 45 arcsec FWHM is displayed 
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Figure 2. Grey-scale of total intensity from the mean of all four L-band images, corrected for attenuation by the primary beam before 
averaging. The resolution is 5.5arcsec FWHM and the grey-scale range, - 1 mjy (beam area) -1 , is marked by the labelled wedge. The 
area is the same as in Fig. ^b). Note that the effects of bandwidth smearing are significant in this image (see Scction l2.2l . 



in Fig. ^b). A higher-resolution (5.5arcsec FWHM) image 
of the same area is shown in Fig. [5] and a detail of the sharp 
bend in the main jet «20 arcmin from the nucleus is plot- 
ted at the same resolution but on a larger scale in Fig. [3] 
Finally, the 5-GHz emission from the inner 4 arcmin of the 
jets is shown in Fig.^Jat a resolution of 2.35 arcsec FWHM. 

We refer to the NW and SE jets as the main and 
counter-jets, as the former is brighter at most distances from 
the nucleus. A striking feature of the main jet is its almost 
constant width between «100 and «400 arcsec from the nu- 
cleus (Fig. Thi s is the collimation shoulder identified by 
IWillis et alJil98ll) : a similar feature is visible in the counter- 
jet, but cannot be traced out as far. The lack of expansion 
over such an extended region is surprising if the jets are con- 



fined solely by thermal plasma associated with the surround- 
ing galaxy group, as a significant pressure gradient would be 
expected on scales of a few hundred arcsec (we argue in Sec- 
tion 15.31 that the core radius of the group-scale plasma is 
~225 arcsec). An alternative possibility is that the jets also 
respond to the J x B forces of their own toroidal fields on 
scales > 100 arcsec. In Section [(J we show that the observed 
polarization structu re is consistent with a dominant toroidal 
component (see also lCanvin et all l2005). but we cannot tell 
from the high-frequency synchrotron emission alone whether 
this component is vector-ordered or has many reversals (ev- 
idence from Faraday rotation is also inconclusive; see Sec- 
tion ^. 31 . The possibility that both pressure confinement and 
magnetic confinement could act together to produce a colli- 
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Table 3. Image resolutions and noise levels. 07 is the off-source 
noise level on the I image; ap the average of the noise levels for 
Q and U . The noise levels were evaluated before correction for 
the primary beam response and apply only at the centre of the 
field for corrected i mages. Asterisks denote the images used by 
ICanvin et afl <2005h 
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Figure 3. Grey-scale of total intensity from the mean of all four 
L-band images, corrected for attenuation by the primary beam 
before averaging and covering the area around the sharp bend in 
the main jet. The resolution is 5.5 arcsec FWHM and the grey- 
scale range is - 0.75 mjy (beam area) -1 . 



matio n shoulder was discussed bv lBridle. Chan fc Henriksenl 
(1981), but they assumed rather different physical parame- 
ters from those we now consider appropriate for NGC 315. 
Our observations show that the collimation shoulder in the 
main jet ends at a bright feature with a sharp edge on 
the side towards the nucleus, inclined by ~60° to the jet 
axis. At this point (marked "Deflection" in Fig. |5J, the flow 
changes direction by ~8° and re-expands with an opening 
angle «10° (defined in terms of t he jet FWHM and the 
angular distance from the nucleus: IWillis et al.l Il98 J) . At 
a similar distance, the counter-jet does not deflect signifi- 



cantly, its surface brightness decreases monotonically away 
fro m the nucleus and it expands less rapidly than the main 
jet iWillis et al.lll98ll Fig.GJ.). 

The brightness distributions in both jets show large- 
scale "banding" - repeated, but irregular alternation of 
bright and faint regions with surface brightnesses differing 
by factors of 1.5 to 2 - along their lengths on arcminute 
scales (Fig. The brightness bands extend across both 
jets but their variations are slower than those in the flaring 
region or at the edges of the jets. These variations could, 
in principle, result either from periods of enhanced activ- 
ity in the nucleus or from interactions between the jets 
and their surroundings. If they were due to fluctuations in 
activity in the nucleus that propagated outwards at con- 
stant velocity /3c, then they would appear at projected dis- 
tances Dj and D c j in the main and counter-jets, respectively, 
where DJD C - = (1 + /3cos0)/(l - /3cos6») and 9 w 38° 
iCanvin et alJ 12005). Any transverse velocity gradients or 
deceleration will complicate this expression and the former 
effect should distort the bands into arcs that are concave to- 
wards the nucleus. We see no obvious relation between the 
distances of the bands in the two jets for any plausible value 
of f3 and no evidence for systematic concave curvature of 
the bands beyond the flaring region. Furthermore, the most 
prominent banding appears to be associated with regions 
where the jets deflect or change their collimation properties. 
It therefore seems more likely that the banding is associ- 
ated with ongoing interactions between the jets and their 
surroundings, although we cannot rule out a contribution 
to large-scale brightness fluctuations from slow variations in 
the jet output. 

The remarkable 180° bend in the main jet at the West 
end of the source is well known from earlier observations. 
Our data (Fig. ^p) show the emission after the bend at a 
resolutio n comparable to th e 610-MHz WSRT image pre- 
sented bv lMack et alJ lll997l) . The brightness distribution at 
the first part of the bend (where the jet deflects by ~100°) 
shows complex structure at 5.5-arcsec resolution. A bright 
ridge runs along the outside edge, with a lane of reduced 
emission next to it (both features are labelled on Fig. 01. 
Note that the sugge stion that the flow is re-energised by an 
intergalactic shock dEnfilin et al.ll200lT) applies to emission 
further downstream, after the second bend. The compact 
knot of emission at the SW edge of the jet just before the 
bend (Fig. |3J coincides in location and shape with the opti- 
cal emission from the flattened background galaxy FGC 0110 
(z = 0.021965) and appears not to be physically related to 
NGC 315. 

The source at RA 00 57 38.710, Dec. 30 22 44.99 (J2000; 
labelled "Background source" in Fig. is unresolved and 
has a flat spectrum. The polarized flux density and E- vector 
position angle vary smoothly across this position, consistent 
with addition of an unpolarized point source to the jet emis- 
sion. There is a faint optical counterpart on the Digital Sky 
Survey and an X-ray point source i s detected at a consis- 
tent p osition with the ROSATPSPC dWorrall fc Birkinshawl 
2000; fig. 2). The source is likely to be a background quasar, 
despite its location on the projected jet axis. There is a 
strong, diffuse component of emission on the axis of the main 
jet ("On-axis enhancement" in Fig. 0^,) at ~225 arcsec from 
the nucleus, with no obvious counterpart in the counter-jet. 

The 2.35-arcsec resolution images illustrate the ini- 
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Figure 4. Grey-scale of total intensity at 5 GHz corrected for attenuation by the primary beam. The resolution is 2.35arcsec FWHM 
and the grey-scale range is - 500 /ijy (beam area) -1 as indicated by the labelled wedge, (a) main jet, (b) counter-jet. 
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tialflaring of the jets (discussed in detail bv ICanvin et alJ 
120051 and Worrall et al., in preparation) followed by recolli- 
mation to a n almost uniform diameter ijWillis et alJll98ll : 
iBridld fl983) Although similar behaviour is observed in 
other sources llLaing fc Bridlel2002al : ICanvin fc Lainel2004) . 
the physical scale on which NGC 315 flares and recollimates 
is unusually large. This is probably a consequence of the low 
external density iWorrall. Birkinshaw fc Hardcastlel [2003 ) . 
and we will explore this idea quantitatively elsewhere, using 
the conservation-law approach developed by Laing & Bridle 
l|2002bl) . The jets in NGC 315 bend slightly as they recolli- 
mate, from a position angle of —48.5° close to the nucleus to 
—52.8° at distances >100 arcsec from the nucleus. The outer 
isophotes of the main and counter-jets are very similar be- 
fore the bend, but the counter-jet is slightly wider at larger 
distances. The main jet is brighter than the counter-jet on- 
axis at all distances from the nucleus, but the counter-jet is 
brighter at the edge between «100 and 200 arcsec. In the flar- 
ing region, by contrast, the main jet is significantly brighter 
than the counter-jet both on-axis and at the edge of the jet 
jCanvin et al]l2005ft . 

We successfully fit the structure of the inner ±70 arcsec 
of the jets in NGC 315 using an intrinsically symmetri- 
cal model in which all apparent differences between the 
main and counter-jets result f rom relativistic aberration and 
beaming dCanvin et alJ 120051) . Such models clearly cannot 
be continued to indefinitely large scales, as FRI jets almost 
always show evidence (e.g. bends and disruption) for asym- 
metric interaction with the external medium. The largest- 
scale structure of NGC 315 shown in Fig. [T[a) is an excel- 
lent example, with jets of unequal length terminating in en- 
tirely different ways. This raises the question: how large is 
the region over which symmetrical, relativistic models can 
be applied? At 400 arcsec from the nucleus, there are clear 
asymmetries in both the deflection and collimation proper- 
ties; these must be intrinsic, although the overall brightness 
asymmetry persists. There is also good evidence for interac- 
tion with the external medium at distances of 70 - 100 arcsec 
on both sides of the nucleus, where both jets bend (symmet- 
rically, this time). It seems likely that intrinsic and relativis- 
tic effects become comparable between 100 and 200 arcsec, 
where the sidedness difference on-axis is the same as on small 
scales, but the edge value is reversed. Asymmetries in ap- 
parent magnetic-field structure, which also imply that the 
flow remains mildly relativistic in this region, are discussed 
in Section HJ Our working hypothesis is therefore that rel- 
ativistic effects dominate the observed differences between 
the two jets only before the first bends, but that environ- 
mental effects become first comparable and then dominant 
at larger distances, although slightly relativistic bulk flow 
probably continues to the largest scales and may remain re- 
sponsible for the generally brighter appearance of the NW 
jet far from the nucleus. 



4 SPECTRA 
4.1 Accuracy 

We define spectral index a in the sense S(u) oc u~ a . We es- 
timated spectral indices both for individual pixels and by in- 
tegration of flux density over well-defined regions. Values at 



5.5-arcsec resolution were determined from weighted power- 
law fits to data at all five frequencies between 1.365 and 
5 GHz. At 1.5-arcsec resolution, spectral indices were calcu- 
lated between 1.413 and 5 GHz. 

There are three main sources of error in the estimate of 
a, as follows. 

(i) The transfer of the amplitude scale from the primary 
calibrator. The errors for the four L-band frequencies are 
likely to be tightly correlated, since they were observed dur- 
ing the same periods, data being taken simultaneously at 
1365 and 1413 MHz, and at 1485 and 1665 MHz. Conse- 
quently, the principal effect of flux-density scale transfer er- 
rors is a constant offset in spectral index. 

(ii) Residual deconvolution effects, typically on scales of 5 
- 20 arcsec. These are approximately proportional to surface 
brightness. 

(iii) Thermal noise. 

We model the error from (ii) as 0.03 times the flux density 
and that from (iii) as the noise level estimated off-source 
(from Table [3J, appropriately integrated. These two contri- 
butions are added in quadrature. In addition, we estimate 
the rms spectral-index offset due to transfer errors in the 
flux-density scale to be 0.02. This should be taken in addi- 
tion to the errors quoted below. 

4.2 Spectral-index images and tomography 

Figs^Ja) and (b) show the spectral indices, a, determined 
from a weighted power-law fit to data at all five frequen- 
cies between 1.365 and 5 GHz at 5.5 arcsec resolution and 
between 1.413 and 5 GHz at a resolution of 1.5 arcsec, re- 
spectively. 

It is clear from Fig.|3a) that there are transverse gradi- 
ents in spectral index where the jets are expanding rapidly. 
These gradients can only be see clearly on spectral-index 
images where the errors are small, and for this reason both 
panels of Fig. are blanked where the rms error in a is 
> 0.05. In order to search for transverse variations over a 
larger area, we need to average along the jets. The gradients 
between 34.5 and 69 arcsec from the nucleus are best dis- 
played by averaging along radii from the nucleus and plot- 
ting the results against angle from the local jet axis, as is 
shown for the main and counter-jets in Fig. [Hfa) . Further 
from the nucleus, where the jets recollimate, we have aver- 
aged along the local jet axis to derive transverse spectral- 
index profiles. The results are shown for two regions in each 
of the main and counter-jets in Fig. |7| The fluctuations in 
these regions are dominated by quasi-periodic deconvolution 
errors: this problem is particularly acute for the counter-jet 
at ~100 arcsec from the nucleus (Fig.^K). The spectral index 
is everywhere consistent with the mean value of (a) = 0.47 
between 70 and 160 arcsec. 

Another method of displaying spatial variations of the 
spectrum is "spectral tomogr aphy" dKatz-Stone fc Rudnickl 
119971: iKatz-Stone et alJll99Sft . This involves the generation 
of a set of images h{r) = I(r, 1/1) — (V2 / vi) at I (r , V2) for 
a range of values of at, where I(r,v) is the brightness at 
position r and frequency v. If the brightness distribution can 
be represented as the sum of two components with different 
spectral indices I(r, v) — S^(v)v~ aa + Sh(r)is~ ab , then the 
"a" component will disappear from the image I t when at = 
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a a . We made a set of images of 7 t with v\ = 1.365 GHz, 
V2 — 5 GHz and a t from 0.4 - 0.7 in steps of 0.01. For 
at ~ 0.44, the outer edges of both jets disappear at distances 
from the nucleus between ~22 arcsec and «80 arcsec. The 
image of It for at = 0.44 is shown in Fig. |H| Further from 
the nucleus, It is close to zero across the whole width of 
both jets for a t ~ 0.47. There is no single value of at for 
which the steep-spectrum component vanishes completely in 
an image of It- 

The main features of the spectral-index distribution are 
as follows. 

(i) Variations in the spectral index are subtle (0.4 < a < 
0.65 everywhere). 

(ii) The spectral index is slightly, but significantly steeper 
in the jet bases than elsewhere. Between 7.5 and 22.5 arcsec 
from the nucleus the mean values at 5.5-arcsec resolution 
are 0.63 and 0.58 in the main and counter-jets, respectively 
(the difference between them is not significant). 

(iii) At 1.5 arcsec resolution (Fig.|SJj), the spectral index 
of the main jet is essentially constant for the first «15 arcsec, 
with a mean spectral index (a) — 0.61, consistent with the 
value determined at lower resolution. 

(iv) Between «15 and 70 arcsec from the nucleus, the 
spectral index is steeper on-axis (a « 0.5) than at the edges 
(a ~ 0.44) in both jets. This is illustrated by the trans- 
verse profiles averaged between 34.5 and 69 arcsec from the 
nucleus (Fig. [SJi) • 

(v) At 1.5-arcsec resolution, the on-axis spectral index is 
slightly higher between 15 and 60 arcsec from the nucleus 
((a) = 0.55; Fig- EJ 3 ) than at smaller distances, a cannot be 
determined to adequate accuracy at the edges of the jet for 
distances >15 arcsec at this resolution. 

(vi) The tomographic analysis shows the spectral gradient 
in a different way: if we subtract off a component with a t = 
0.44 (Fig. [HJl, the emission at the edges of the jet and at 
large distances from the nucleus essentially vanishes. What 
remains (positive in Fig. |HJ corresponds to the jet bases and 
to a ridge of steeper-spectrum emission at larger distances. 
The latter is clearly visible in both jets. 

(vii) The flatter-spectrum edge first becomes detectable 
at ~15 arcsec from the nucleus and widens thereafter. It oc- 
cupies the entire width of the jet from «70 arcsec outwards. 
The transition between steeper and flatter spectrum on-axis 
is poorly defined. 

(viii) There is no evidence for any transverse spectral gra- 
dient at larger distances, after the jets recollimate, although 
the data are noisy and do not cover quite the full width of 
the jets (Fig. 

(ix) This is confirmed by the tomographic analysis: It for 
the outer parts of the region vanishes for a t = 0.47, the 
mean spectral index, confirming that a is constant within 
our errors. 
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Figure 5. (a) False-colour plots of spectral index, a over the 
range 0.3 - 0.8. (a) a determined by weighted power-law fits to 
data at all five frequencies between 1.365 and 5 GHz. The res- 
olution is 5.5 arcsec FWHM. Data are plotted only where the 
rms error in a (excluding flux-density scale offsets) is <0.05 and 
the signal-to-noise ratio at all frequencies is >5. All of the im- 
ages were made using a zero-spacing flux density and the (small) 
residual zero-point offsets were subtracted before fitting. The 
75 fi Jy (beam area) _ 1 contour from the mean L-band / image 
is also plotted, (b) a between 1.413 and 5 GHz at a resolution 
of 1.5 arcsec FWHM for the area shown by the box in panel (a). 
Data are plotted only where the rms error in a (again exclud- 
ing flux-density scale offsets) is <0.05 and the difference between 
estimates of a before and after zero-level subtraction <0.01. Rep- 
resentative contours of total intensity at 5 GHz are also plotted. 



4.3 Deprojection of the spectral-index 
distribution 

iKatz-Stone et alJ l|l999ll and lKatz-Stone fc Rudnickl lll997D 

suggested that the spectral index of an on-axis component 
in a jet is the value of a t at which the component appears 
to vanish against the background of the surrounding emis- 
sion (exactly as for an edge component such as that in the 



NGC 315 jets; Section r4.21 and can therefore be derived sim- 
ply from a tomographic analysis. This requires an additional 
assumption which may not be correct, namely that when the 
the on-axis component is subtracted, the remaining emission 
has a smooth brightness distribution (this is not true for the 
model proposed below). For NGC 315, our th ree-dimensional 
model of the emissivity iCanvin et al.ll2005T) gives a good fit 
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Figure 6. (a) Averaged transverse profile of the spectral in- 
dex determined from a weighted power-law fit to data at all 
five frequencies between 1.365 and 5 GHz, as in Fig. EJa). The 
resolution is 5.5 arcsec FWHM. The data were averaged along 
radii from the nucleus between 34.5 and 69 arcsec and are plot- 
ted against angle from the jet axis, here taken to be in posi- 
tion angle —48.5°. Full curve: main jet, dashed curve: countcr- 
jet. (b) Profil e of jet /counter-jet sidedness ratio, from the data in 
Canvin et alJ l2005fl , averaged as in panel (a) , but with a resolu- 
tion of 2.35 arcsec FWHM. 

to the observed emission at 5 GHz, so we isolated the on-axis 
component and measured its spectral index, as follows. 

(i) We first used the tomography image with a t = 0.44 
as a template for the on-axis component. To a reasonable 
approximation, this defines a cone with a half-opening angle 
of 13.6° projected on the sky. Assuming a n angle to the line 
of sight of 9 = 38° JCanvin et al]l200Sft . the half-opening 
angle in the jet frame is 8.3° . 

(ii) We then m ade a model 5-GHz image, as in 
ICanvin et al.l i2005l) but with the emissivity set to zero 
within this cone, convolved it to a resolution of 5.5 arcsec 
and subtracted it from the observed 5-GHz / image. The 



Figure 7. Averaged transverse profiles of the spectral index de- 
termined from weighted power-law fits to data at all five frequen- 
cies between 1.365 and 5 GHz, as in Fig. |3a). The resolution is 
5.5 arcsec FWHM. The data were averaged along the local jet 
axis, taken to be in position angle —52.8°. Full curve: main jet; 
dashed curve: counter-jet. The horizontal lines indicate the mean 
spectral indices for the regions, which are close to 0.47 in all cases, 
(a) 69 - 112.5 arcsec, determined by averaging the spectral-index 
image, (b) 114 - 157.5 arcsec, determined by summing individual 
/ images and making weighted power-law fits (the spectral-index 
image is too noisy to be averaged towards the jet edges). 

residual corresponds to the observed emission from within 
the central cone. 

(iii) We scaled the model to 1.365 GHz assuming a spec- 
tral index of 0.44, as appropriate for the edge emission and 
then subtracted it from the observed 1.365-GHz image. 

(iv) Both residual images showed little emission towards 
the edges of the jets, implying that the model subtraction 
was reasonably accurate. 

(v) Finally, we derived a spectral-index image for the on- 
axis component alone. This is shown in Fig. El 

Fig. El shows that the spectrum of the on-axis compo- 
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Figure 8. The grey-scale plot shows a difference image — 
{v 2 /vi) at I{v 2 ) for u\ = 1.413 GHz, v 2 = 5 GHz and at = 0.44 
over the range —1 to +lmjy(beam area) -1 . The resolution is 
5.5arcsec FWHM. A few contours of the 5 GHz / image are su- 
perposed to outline the jet structure. 
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Figure 10. A sketch of our proposed three-dimensional distri- 
bution of spectral index in the jets of NGC 315. The sketch is 
in a plane contai ning the jet axis, as suming an angle to the line 
of sight of 37.9° JCanvin et al.ll2005l) and a linear scale is given. 
10 kpc along the jet projects to 18arcsec on the plane of the sky. 
The values of a at various points in the jet are marked, together 
with the approximate extent of the steep-spectrum component 
(dott ed) and the region s of high emissivity and rapid decelera- 
tion iGanvin et al.ll2005l) . 



they are slightly higher than the values measured on-axis 
before subtraction, which include a contribution from the 
flatter-spectrum, off-axis component. We conclude that the 
structure observed in Fig. |SJa) does not result entirely from 
the superposition of two components with constant, but dif- 
ferent spectral indices. The sketch in Fig. llOl summarizes our 
results on the distribution of spectral index in the jets. 
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Figure 9. The false-colour plot shows an estimate of the spectral 
index for the on-axis emission component alone in the range 0.4 Sj 
a $C 0.8. The spectral-index distribution is truncated at a distance 
of 66.5 arcsec along the axis and is shown only where the residual 
surface brightnesses exceed 0.5 and 1 mjy (beam area) -1 at 5 and 
1.365 GHz, respectively. A few contours of the 5 GHz / image arc 
superposed to outline the jet structure. The area covered is the 
same as that in Fig. Inland the resolution is 5.5 arcsec FWHM. 

nent flattens slightly with distance from the nucleus. The 
mean values of a between 7.5 and 22.5 arcsec are 0.60 and 
0.61 for the main and counter-jets, respectively. The cor- 
responding values for distances between 22.5 and 66.5 arc- 
sec are 0.56 and 0.55. For the main jet, these values are 
consistent with the measurements at 1.5-arcsec resolution 
without subtraction. At 5.5-arcsec resolution in both jets, 



4.4 Comparison with other sources 

Typical spectral indices for the flarin g regions of other FRI 
radio jets are in the range 0.5 - 0.6 llYoung et al. 2005 and 
references therein). Indeed, lYoung etHuTT^OSl) suggested 
that FR I jets have a "canonical" low-frequency spectral in- 
dex of a = 0.55, but this conclusion is based primarily on 
lower-resolution data than we consider here. Their canoni- 
cal value is intermediate between that for the jet bases in 
NGC 315 and the significantly flatter spectra seen at larger 
distances. 

There are very few published studies of spectral vari- 
ations in the flaring regions of FRI jets. Short regions 
of slightly steeper spectrum than the average have been 
detected in the bases of t hree other FR I jets: 3C 449 
jKatz-Stone fc Rudnickll997l) . PKS 1333-33 jKilleen et al. I 
1986) and 3C 31 (Laing et al., in preparation). The measure- 
ment of a = 0.7 ± 0.2 in 3C449 applies to the faint inner 
jets (corresponding to the innermost 5 arcsec in NGC 315), 
but differs marginally from that of the brighter nearby emis- 
sion, given the large error. In PKS 1333—33, the spectrum 
flattens from a ~ 1.0 to a ~ 0.6 between 10 and 35 arcsec 
from the nucleus (2.4 - 8.4 kpc in projection); this includes 
both the faint inner jets and the bright part of the flaring 
region, as in NGC 315. In 3C 31, there is a steeper-spectrum 
region extending to «6 arcsec in the main jet (plausibly also 
in the counter-jet). 

There is a flatter-spectrum edge on one side of the main 
jet in 3C31 (Laing et al., in preparation). Transverse vari- 
ations of spectral index in the sense that the spectrum is 
steeper on-axis have not been reported in any other sources, 
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but their flaring regions are smaller in both linear and an- 
gular size and are difficult to resolve. Gradients at the level 
of Aq ~ 0.05 - 0.1 are also tricky to detect without data at 
more than two frequencies. 

The tendency for the jet spectrum to be flatter at the 
edge is in the opposite sense to the spectr al gradients found 
on large scales in tailed radio so urces llKatz-Stone et all 
Il999t iKatz-Stone fc Rudnicldll997t Laing et al., in prepa- 
ration) but the latter effect occurs in completely different 
regions of the jets, where they merge into the tails. 

4.5 Acceleration mechanisms 

The X-ray emission de tected by 

IWorrall. Birkinshaw fc Hardcastld J2003D in NGC315 
coincides with the steeper-spectrum (a = 0.61) region at 
the base of the main jet. The form of the synchrotron 
spectrum in those FR I jet bases which emit X-rays i s now 
well established (Qj ardcastlc, Birkinshaw &^Worral|_ 2001; 
|Hardcastle_ej i _ai|]2j)0^ alJ l2003t lHardcastle et al.l 

1200.4 IPerlman fc Wilson! l2005h . It can be characterised 
as a broken power law with spectral indices of 0.5 - 
0.6 at radio wavelengths and 1.2 - 1.6 in X-rays. The 
spectrum of NGC 315 is consistent with this pattern, 
although its high-frequency slope is po orly constrained 
iWorrall. Birkinshaw fc Hardcastlel 120031) . This spectral 
shape is consistent with synchrotron emission from a 
single electron population; ongoing particle acceleration is 
therefore required (see also Worrall et al., in preparation). 
The magnetic-field st rengths calcu la ted f or the on-axis 
emissivity model of ICanvin et al] J2005) , assuming a 
minimum-pressure condition, range from 3.3 nT at 6arcsec 
in projection to 0.44 nT at 69arcsec. For electrons with 
Lorentz factor 7 radiating at the synchrotron critical fre- 
quency v c in a magnet ic field B, we have f c /Hz = 41.99 7 2 
(B/nT) llLongairll994h . so at 5 GHz, 6 x 10 3 < 7 < 1.6 x 10 4 
and in the X-ray band at 1 keV, 4 x 10 7 < 7 < 1.1 x 10 s . 

The flatter-spectrum edges occur where ICanvin et al] 

(2005) infer substantial velocity shear across the jets. If the 
jets are relativistic and faster on-axis than at their edges, 
the approaching jet always appears more centre-brightened 
than the receding one. This difference is reflected in the 
sidedness-ratio image. The average transverse sidedness pro- 
file between 34.5 and 69arcsec from the nucleus is shown in 
Fig. HJb) for comparison with the spectral-index profile over 
the same region. The velocity profile is modelled as a trun- 
cated Gaussian fu nction with 3 = 0.3 8 on-axis and 0.22 at 
the edge, although lCanvin et al] i2005T) note that the on-axis 
velocity may be larger (0 « 0.5) if the shear occurs over a 
narrow range of radii in the jet: this might give a better fit 
to the sidedness profile. 

The flatter-spectrum edge first becomes visible 
«15 arcsec from the nucleus. This coincides to within the 
errors with: 

(i) the start of rapid deceleration, as inferred by 
ICanvin et al] (|200jf), 14arcsec in projection from the nu- 
cleus; 

(ii) the end of the region of enhanced radio and (in the 
main jet only) X-ray emissivity dCanvin et al]l2005l Worrall 
et al., in preparation). 

(iii) The first point at which the observed jet/counter-jet 



sidedness image gives any evidence f or transverse veloc ity 
gradients, 16arcsec from the nucleus llCanvin et al.l |2005). 

The deceleration and enhanced emission regions are marked 
on Fie. 1101 Note that the detection of transverse gradients in 
sidedness and spectral index may be limited by resolution. 

The changes of spectral index observed in NGC 315 sug- 
gest that (at least) two different electron acceleration mech- 
anisms are required, as follows. 

(i) The first mechanism dominates at the base of the flar- 
ing region (the initial 15arcsec in NGC 315) and may con- 
tinue at a lower level on the axis of the jet to «70 arcsec. 
It generates emission from radio to X-ray wavelengths and 
has a characteristic spectral index a w 0.6 in the former 
band. Three pieces of evidence suggest that this mechanism 
is dominant where the jet is fast (0 ss 0.9). Firstly, its char- 
acteristic spectral index is observed across the whole of the 
jet width in NGC 315 u ntil the start of rapid deceleration. 
Secondly, in both 3C 31 jLaing fc Bridlel2004l) and NGC 315 
(Worrall, Birkinshaw & Hardcastlc 2003), the bright X-ray 
emission occurs upstre am of the deceleration region. Finally, 
lLaing fc Bridld i2004l) show from radio data alone that sig- 
nificant injection of fresh relativistic particles is required be- 
fore the start of deceleration in 3C 31 to counter-balance 
adiabatic losses. 

(ii) The second mechanism causes the flattening of the 
spectrum towards the edges of the jets observed from 
15 arcsec outwards, but eventually spreading over the en- 
tire jet width. It produces spectral indices in the range 
0.44 < a < 0.5 for electrons emitting at radio wave- 
lengths in NGC 315 and appears to be associated with ve- 
locity shear across the jets. A possible candidate for the this 
mechanism is the shear acc eleration process described by 
iRieger fc Duffvl (|2004 12005I) . but their estimates of the ac- 
celeration timescale for electrons in conditions appropriate 
to FR I jets are very long (at least if the mean free path ~ 
gyro-radius) , and it is unclear whether the process is efficient 
enough to influence the spectrum. 



5 FARADAY ROTATION AND 
DEPOLARIZATION 

5.1 Faraday rotation 

In order to investigate the variations of Faraday rotation 
along the jets of NGC 315, we made images of rotation mea- 
sure (RM) at a resolution of 5.5 arcsec by least-squares fit- 
ting to the relation x(^ 2 ) = X(°) + RMA 2 (where x is the E- 
vector position angle) for all 5 frequencies between 1.365 and 
5 GHz. The fits were weighted by errors in \ derived from 
Table El We excluded a small region around the core which 
was affected by residual instrumental polarization and in- 
cluded only points where the rms error in position angle was 
<15° at all frequencies. The resulting RM image is shown 
in Fig. EI a). The fit to a A 2 law is very good everywhere: 
two examples are shown in Fig. 1121 The extreme value of 
RM » — 90radm -2 results in rotations of » 1° between the 
two centre frequencies combined in our 5 GHz dataset and 
«0.4° and ~5° across the bands at 5 and 1.4 GHz, respec- 
tively. The worst of these effects, rotation across the band 
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Figure 11. False colour images of rotation measure for the jets in NGC 315. In panels (a)- (c), the RM is derived from fits to E-vector 
position angles at 5 frequencies, (a) Observed RM, in the range -90 to -60 radm -2 . The crosses mark the positions corresponding to 
the x ~ ^ 2 plots in Fig. 1 1 21 (b) RM after subtraction of a linear function of distance along the jet derived from an unweighted fit to 
the data in the range —15 to +15 radm" 2 . (c) As (b), but for a fit to gradients along and across the jets, (d) RM at a resolution of 
1.5arcsec FWHM, derived from the position-angle difference between 1.413 and 5 GHz, using the lower-resolution RM image to resolve 
mr ambiguities. Data are plotted only where the rms error in RM <2radm -2 . The area covered by this figure is indicated by the box 
on panel (a). 



at the lowest frequencies, results in a spurious depolariza- 
tion <0.1%, which is negligible compared with errors due 
to noise. The images of Q and U show little power on large 
spatial scales and our estimates of position angle should be 
reliable over the full range of distances from the nucleus 
shown in Fig. Illf a'l. The area over which the RM can be 
determined accurately is limited primarily by the primary 
beam at 5 GHz (540arcsec FWHM). 

The mean RM is — 75.7radm -2 . The most obvious fea- 
ture of the RM image is a nearly linear gradient along the 
jets, as shown by the profile in Fig. I13f a1 . In order to reveal 
smaller-scale structure in the RM, we initially fitted and sub- 
tracted a function RM = RMo + a x x, where x is measured 
along the axis from the nucleus (a x and RMo are constants). 
We used an unweighted least-squares fit, as any attempt to 
weight by the estimated errors caused the brightest part of 



the main jet to be fitted at the expense of other regions. 
The gradient is a x — 0.025 radm -2 arcsec -1 and the resid- 
ual image is shown in Fig. I13f b~l . This indicates that vari- 
ations in RM across the jet are also significant, as can be 
seen more clearly in averaged profiles, particularly between 
69 and 112.5 arcsec from the nucleus (Fig. H4( d'). Closer to 
the nucleus, the gradient is barely visible (Fig. 114b ). but the 
jets are narrower there and the profiles are consistent with 
the gradient measured at larger distances. The transverse 
variation also appears to be linear, so we fitted and sub- 
tracted a function RM = RMo + a x x + a y y, where y is a 
coordinate transverse to the jet and a y is a constant. The 
gradients along and transverse to the jet axis become a x = 
0.018 radm -2 arcsec -1 and a y = 0.051 radm -2 arcsec -1 , re- 
spectively. Taken at face value, the best-fitting gradient is 
0.054 rad m -2 arcsec -1 at an angle of 72° to the mean jet 
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Figure 12. Plots of E-vector position angle, \, against A 2 for 
representative positions in the main jet (panel a) and counter-jet 
(panel b). The coordinates are given on the plots and the positions 
are those marked by crosses in Fig. Illl a). 



axis. Note, however, that the transverse gradient is essen- 
tially determined by a subset of the data from the widest 
parts of the main and counter-jets (Fig. I14b 1. 

Removal of the large-scale gradient leaves fluctuations 
in the local mean RM which appear significantly larger on 
the counter-jet side (Fig. Illb ). By definition, the signal-to- 
noise ratio in I is lower in the counter-jet. Although this is 
partially offset by a higher average degree of polarization, 
the errors in RM are still larger than in the main jet. To 
evaluate the significance of the fluctuations, we considered 
only points with fitting errors ^2.5 rad m -2 and calculated 
the expected errors in the means for boxes of length 30 arcsec 
along the jet axis containing more than 50 such points (the 
errors are corrected for oversampling) . The results are shown 
in Fig. I13f b1. The fluctuations are significant, and form an 
ordered pattern with a typical scale ~100 arcsec. They are 
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Figure 13. Variation of the mean and rms RM along the jets. 
The average jet axis is taken to be at position angle —52.8°. (a) 
A profile of RM averaged over boxes of length 30 arcsec along the 
jet axis, (b) Mean RM profile, after subtraction of a linear fit. 
Only points with fitting error ^2.5 rad m — 2 were used and boxes 
are only plotted if they contain more than 50 such points. The 
vertical bars indicate the error on the mean, (c) the rms RM, 
with respect to the mean of the points in the box. Dashed bars 
denote the uncorrected values, URMraw ; full bars show the values 
after making a first-order correction for fitting error, <r RM (see 
text). The curve plotted in panels (b) and (c) is the model rms 
RM profile described in the text. 
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Figure 14. Transverse variation of RM after subtracting a uni- 
form gradient along the axis for the main jet (full line) and 
counter-jet (dashes), (a) Averages along radii from the nucleus 
between 34.5 and 69 arcsec plotted against the angle from the 
axis, which is taken to be in position angle —48.5°. (b) Averages 
along the jet (in position angle —52.8°) between 69 and 112.5 
arcsec, plotted against distance from the axis. 



larger by a factor of ~2 in the counter-jet and the first bin 
of the main jet (within 30 arcsec of the nucleus) compared 
with the rest of the main jet. 

The residual fluctuations within the boxes (i.e. after 
subtracting the local mean) are comparable with the errors 
in RM except in the brightest regions close to the nucleus. 
We therefore made a first-order correction to the rms RM, 
ORMraw, by subtracting the fitting error ant m quadrature 
to give g rm = (o-RMraw ~ it) 1/2 - The profiles of urm™ 
and ctrm are both plotted in Fig. I13f c). for the same se- 
lection of points as in Fig. Il^f bl . The corrected profile is 
very uncertain, but suggests that orm has a maximum of 
2radm -2 close to the nucleus on the counter-jet side and 
may be slightly asymmetric in the sense that the rms RM 
is lower in the main jet than the counter-jet at the same 
distance from the nucleus. 

We can image these smaller-scale fluctuations directly 
only at the bright base of the main jet. There, the RM 
at a resolution of 1.5 arcsec FWHM can be derived accu- 



rately from the difference between 1.413 and 5 GHz position- 
angle images, using the lower-resolution data to resolve the 
nit ambiguity. Fig. Illf d) shows the RM at this resolution. 
Data are plotted only where the rms error in the fitted RM 
<2radm~ 2 . Fluctuations are clearly detected: the rms is 
(TRMraw = 2.1radm -2 , giving ctrm = 1.6radm -2 after mak- 
ing a first-order correction for fitting error, as above. This 
is in good agreement with the value for the innermost bin 
of the profile of rms RM for the main jet at lower resolution 

(Fig. Est). 



5.2 Depolarization 

The variation of p with wavelength at low resolution po- 
tentially measures fluctuations of RM across the observing 
beam which cannot be imaged directly with adequate sensi- 
tivity. This variation is small, and is best quantified by fit- 
ting to the first-order approximation p(A 2 ) « p(0) ±p'(0)A 2 
where p'(A 2 ) = dp/d(\ 2 ). We expect p'(0) < (depolar- 
ization) under most circumstances. The quantity p'(0)/p(0) 
is directly related to the commonly quoted depolarization 
ratio, but is biased in the present case, since both the gradi- 
ent and the degree of polarization depend directly on a sin- 
gle high-frequency measurement (at 5 GHz), so deviations in 
p(0) and p'(0) are anticorrelated. We therefore tested for the 
presence of depolarization using the gradient p'(0) alone. We 
derived p'(0) by weighted least-squares fitting to images of 
p at the 5 frequencies between 1.365 and 5 GHz at a resolu- 
tion of 5.5 arcsec FWHM. Two sets of p images were used: in 
the first, an estimate of the local zero- level was subtracted 
from the I images before calculating p = P/I; in the sec- 
ond, the original I images were used (any differences indicate 
systematic errors in the estimation of large-scale structure) . 
The fitting weights were the inverse squares of errors in p 
derived from the values in Table |H] and points were only 
included if the errors were <0.3 at all frequencies. The re- 
sulting polarization gradients are very small, and there are 
no obvious variations. In order to determine the significance 
of the gradients, we measured their mean values over the 
main and counter-jets. The maximum scale of structure im- 
aged accurately in total intensity is «300 arcsec (Table|2J), so 
we calculated the means for points between 9 and 150 arcsec 
from the nucleus along the jet axis (excluding the small re- 
gion around the core to avoid spurious instrumental polar- 
ization, as in Section 15.11 . The mean values for the main 
and counter-jets derived from the zero-level corrected im- 
ages were <p'(0)) = 0.03 ±0.08 and 0.16 ±0.16, respectively. 
Without zero-level correction, the values become 0.05 ±0.08 
and 0.22 ± 0.16. We conclude that the increase of p with 
wavelength for A ^ 0.22 m {v 1.365 GHz), which is in any 
case opposite to the expected effects of Faraday rotation, is 
not significant. 



5.3 The origin of the rotation measure 

NGC315 has Galactic coordinates I = 124.°6, b = 
—32.° 5. This is on the ou t skirts of Region A of 
ISimard-Normandin fc Kronberj {l980), where the major- 
ity of sources have large nega tive RM's (~ — l OOrad m -2 
at the centre of the region). iDineen fc Colesl <|2005T) de- 
rived spherical-harmonic models of the Galactic RM distri- 
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bution by fitting to RM's of large numbers of extragalac- 
tic sources. These models predict a Galactic contribution 
?s -47radm~ 2 at the position of NGC315. The bulk of the 
mean RM of — 75.7radm -2 is therefore likely to be Galactic 
in o rigin. 

ISimonetti fc Cordesl lll98rJ) determined RM variations, 
which they argued to be primarily Galactic, across a num- 
ber of sources within Region A (their Region 1). Their plot 
of squared RM difference ARM 2 against separation suggests 
<ARM 2 } 1/2 ~ 10 - 30radm~ 2 on a scale of 700arcsec, but 
with large uncertainties. The observed ARM (« 15radm -2 
on the same scale along the jets of NGC 315; Fig. I13b ^l 
is again consistent with a Galactic origin. The maximum 
change in RM across the jet (~ 4radm -2 over 80arcsec; 
Fig. 114b ) is within the range o f the upper limits plot- 
ted bv lSimonetti fc Cordesl dl986h . We note, however, that 
NGC 315 is located outside the core of Region A, so it may 
be that smaller values of (ARM 2 ) 1 / 2 are appropriate. 

The linearity of the RM gradient along the jet and the 
fact that the maximum gradient is aligned neither with the 
jet axis nor with the minor axis of the galaxy both imply that 
little of the Faraday rotating medium is associated with the 
jets or with the host galaxy. In either case we would expect 
some non-linear variation with distance from the nucleus. 

We conclude that most of the mean RM and its linear 
gradient are likely to be Galactic in origin, but a significant 
contribution from material local to NGC 315 is not ruled out. 
In particular, we cannot exclude the hypothesis that some 
of the RM gradient transverse to the jet results from an or- 
dered toro idal field within or just outside the jet IlLaing 1981; 
lBlandfordlll993T) . The associated position angle rotation be- 
tween 1.365 and 5 GHz is at most a few degrees and signif- 
icant departures from A 2 rotation would not be detectable 
even if the thermal plasma responsible for the Faraday rota- 
tion is mixed with the synchrotron-emitting material jBurnl 
1966). If local toroidal fields are solely responsible for the 
transverse RM gradients, then their vector directions must 
be the same in the main and counter-jets. 

The residual RM fluctuations are qualitatively very sim- 
ilar to those in 3C31 (Laing et al., in preparation) but have 
amplitudes that are 10 times smaller on similar angular 
scales. The larger-scale (~100 arcsec) fluctuations are sys- 
tematically lower on the main (approaching) jet side. The 
distribution of fluctuations on smaller scales is also consis- 
tent with such an asymmetry, but is not well determined. 
As with the transverse gradients discussed earlier, the ob- 
served posit ion- angle rotations are too small to be sure that 
the RM fluctuations are due to foreground plasma, but the 
asymmetry between approaching and receding jets suggests 
an origin in a distributed magnetoionic medium surrounding 
the host galaxy, a possibility we now explore. 

The thermal plasma associated with NGC 315 
and observed using Chandra can be described by 
a beta model with a core r adius of 1.55 arcsec 
iWorrall. Birkinshaw fc Hardcastlel I2003T) . This cannot 
be responsible for the RM fluctuations, which occur on far 
larger scales. The most likely hot plasma component to be 
responsible for the Faraday rotation would be associated 
with the poor group of galaxies surrounding NGC 315 
jNoltheniuslll993l : iMiller et al]l2002l) . but has not yet been 
detected in X-ray observations. Since RM fluctuations are 
seen in both jets, a spherical distribution is plausible. In 



order to make a rough estimate of the parameters of the 
putative group component, we took a simple model for the 
field structure in which cells of fixed si ze I at radius r con- 
tain r andomly orie ntated fields B(r) llLawler fc Dennisonl 
Il982t iFeltenl Il99rj . The density distribution was taken 
to be a beta model: n(r) = n (l + r 2 / r 2)-3/3atm/2 wit]l 
B(r) oc n(r) N . We derived the variation of Urm along the 
projection of the jet axis by numerical integration, assuming 
that the jets have 9 = 37.9° everywhere. This is a similar 
approach to the calculation o f dep olarization asy mmetry by 
iGarrington fc Conwavl dl99ll) and lTribblei dl992[) ; our code 
also reproduces the analytical results ofjFelten|_ ll 19961) for 
6 = 9 ° and N = or 0.5. Following iDolag et alj <200ll) : 
iDolad l)2006l) . we assumed that B(r) oc n(r) or N = 1 
(an empirical result derived from RM's of radio sources in 
and behind cluster and rich groups). Our detection of RM 
variations on a range of scales is qualitatively consistent 
with the idea that the power spectrum of the magnetic-field 
fluctuations is a power law iTribbld Il99lt lEnfilin fc Vogd 
l2003tlMurgia et al .120041 : 1 Vogt fc Enfilinl2005h . but our data 
are to o noisy and poorly sa mpled to constrain its functional 
form. iMurgia et alJ 1120041) show that a single-scale model 
gives a very similar relation between the RM variance 
and radius r to that derived for the more realistic case of 
a power-law power spectrum provided that I is interpreted 
as the correlation length of the magnetic field. Finally, 
we fitted the resulting orm curves by eye to the profiles 
of RM fluctuations on different scales in Fig. I13f b) and 
(c). We fixed the value of /3 a tm = 0.5 and adjusted the 
core radius to give a reasonable fit to the profile. For both 
plots, we found that r c m 225 arcsec gave an adequate fit. 
Strictly speaking, Orm is the RM variance evaluated over a 
window much larger than the maximum fluctuation scale, 
whereas the profiles in Fig. 1131 describe fluctuations over 
two different ranges of scale. We have therefore added the 
two model variances to give a rough estimate of the total 
(the curves shown in Fig. 113b and c actually have the same 
amplitude). Note that the process of removing a linear 
trend from the RM profile will have suppressed some power 
in large-scale fluctuations, particularly transverse to the 
jet. The amplitude of the model variance profile is related 
to the central dens ity and field, and to the correlation 
length lFeltedll996h : (n /m- 3 ) 2 (B /nT) 2 (//kpc) « 700. 
This is a very rough estimate, but is enough to establish 
that a low-density group-scale gas component with a low 
magnetic field can generate the observed Faraday rotation. 
Plausible parameters mig ht be n ~ 600 m" 3 , B w 
0.015 nT (0.15uGauss) and I ^ lOkpc. We note that the 
recollimation of the jets may also require a large-scale hot 
gas component. 



6 MAGNETIC-FIELD STRUCTURE 

Vectors whose magnitudes are proportional to p at zero 
wavelength and whose directions are those of the apparent 
magnetic field inferred from the rotation-measure fit of Sec- 
tion l5.1l are plotted in Fig. 1151 At higher resolution, we de- 
rived the apparent field direction by interpolating the RM 
image onto a finer grid and using it to correct the observed 5- 
GHz position angles fFig. I16H . The apparent field structure 
in the flaring region (up to ~70 arcsec from the nucleus) 
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Figure 15. Vectors whose lengths are proportional to the degree of polarization at zero wavelength, p(0), from the analysis of Section l5,2l 
and whose directions are those of the apparent magnetic field, derived from the 5- frequency RM fit of Section 15, II The vector scale is 
indicated by the labelled bar. The grey-scale is of total intensity at 4.985 GHz and the resolution is 5.5arcsec FWHM. 



is discussed extensively by ICanvin et all l)200Rl) . Here, we 
concentrate on larger scales, after the jets recollimate. 

The signal-to-noise ratio for individual points, partic- 
ularly at the jet edges, is often <3 in linear polarization, 
causing them to be blanked in Figs ll5l and ll6l In particular, 
it is impossible to see whether the parallel-field edge contin- 
ues to large distances in the main jet. We therefore adopted 
the following procedure to derive the average degree of po- 
larization. 

(i) We first corrected the observed 5-GHz, 2.35-arcsec po- 
sition angles for Faraday rotation using the linear model de- 
rived in Section lb.il which is defined everywhere in the field, 
unlike the RM image. 

(ii) We then changed the origin of position angle to be 
along the jet axis, so that apparent field along or orthogonal 
to the jet appears entirely in the Q Stokes parameter (also 
verifying that there is very little signal in U). 

(iii) We then integrated Q and I along the jet axis from 
69 - 113 and from 113.5 - 157.5 arcsec from the nucleus (the 
same areas used for the profiles of spectral index in Fig.0. 
Lack of short spacings at 5 GHz precludes extension of this 
analysis to larger distances. 



(iv) Finally, we divided the results to give transverse pro- 
files of Q/I. Provided that the apparent field is either along 
or orthogonal to the jet axis, p = \Q/I\. We have chosen the 
sign convention so that Q > for transverse apparent field 
and < for longitudinal field. 

The resulting profiles are shown in Fig. 1171 Figs 1151 - 1171 
show that the apparent field configuration found in the flar- 
ing region - transverse on-axis and longitudinal at the edges 
- persists until at least 160 arcsec in both jets. In particular, 
the longitudinal-field edge of the main jet is easily detected 
in the profiles, even though it is not clearly visible on the 
images. The main difference fro m the corresponding profiles 
for the flaring region (fig. 9 of ICanvin et al.ll2005l) is that 
the on-axis polarization is higher in both jets at larger dis- 
tances. This is a continuation o f the tr end in the longitudinal 
profile shown bv lCanvin et alJ <2005l their fig. 8). As in the 
flaring region, the on-axis (perpendicular) polarization is al- 
ways higher in the counter-jet, reaching levels close to the 
theoretical maximum of po = 0.69 for the observed spectral 
index. In the main jet, p « 0.4 on-axis. Both jets are very 
hig hly polarized a t their edges. 

ICanvin et alJ ll200fJl modelled the three-dimensional 
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structure of the field in the outer parts of the flaring re- 
gion as a mixture of toroidal and longitudinal components of 
roughly equal magnitude on-axis but with the former domi- 
nant at the edge of the jet. The differences between the two 
jets are attributed to relativistic aberration, so the fact that 
they persist after the jets recollimate suggests that there 
is little further deceleration in this region, despite the re- 
versal in sidedness at the edges of the jets (Section [3J . A 
decrease in the on-axis longitudinal field component, bring- 
ing the configuration closer to a purely toroidal one, would 
result in a polarization distribution consistent with that ob- 
served. Some longitudinal component must remain, however, 
otherwise p would be close to po on the axis for both jets. 
The very high degree of polarization at the jet edge requires 
the radial field com ponent to be very sm all, as inferred for 
the flaring region bv lCanvin et alJ i2005l) . 
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7 SUMMARY 

We have imaged the jets in the nearby FRI radio galaxy 
NGC 3f 5 with the VLA at five frequencies in the range 
f .365 - 5 GHz and at resolutions ranging from 45 - f .5 arcsec 
FWHM. Our total intensity observations reveal new details 
of the structure, particularly around the sharp bend in the 
main jet. 

The flaring regions of both jets, where they initially 
expand rapidly and then recollimate, show a complex and 
previously unknown spectral structure. Within 15 arcsec of 
the nucleus, the spectral index has a uniform value of a = 
0.61 in both jets. This region is associated with strong X-ray 
emission in the main jet, high radio emissivity, complex fila- 
mentary structure, and fast flow with 8 ~ 0.9 iCanvin et al] 
l2005tlWorrall. Birkinshaw fc Hardcastlell2003lh Between 15 
and 70 arcsec, the spectrum is steeper on-axis than at the 
edges of the jet. We have developed a novel deprojection 
technique which allows us to isolate two spectral compo- 
nents. The first (on-axis) forms a continuation of the jet 
base, its spectral index flattening gradually from 0.61 to 
0.55. The second (at the edge of the jet) has a w 0.44 and 
is associated with a region where strong shear is inferred 
JOanvin et aljEoOSh . We speculate that two different accel- 
eration mechanisms are involved, one associated with fast 
flow, dominant close to the nucleus and capable of acceler- 
ating electrons to the very high energies required to produce 
X-ray emission (7 ~ 10 8 ), the other being driven by shear 
and generating the flatter spectral indices seen at the edges 
of the jet. Both mechanisms must efficiently generate the 
electrons with 7 ~ 10 4 which radiate at cm wavelengths. At 
distances >70 arcsec, the spectral index is consistent with 
a « 0.47 everywhere. 

We have imaged the variations of Faraday rotation over 
the jets. All of the rotation is resolved and must originate 
mostly in foreground material. There is no detectable de- 
polarization. The largest contributions - a constant term 
and a linear gradient - are probably Galactic in origin. We 
have also detected residual fluctuations of «1 - 2 radm -2 
rms on scales ~ 5 - 100 arcsec. The amplitude of fluctua- 
tions on scales > 30 arcsec is larger by a factor ~2 for the 
counter-jet, consistent with an origin in magnetoionic mate- 
rial around the source, but not in the known X-ray-emitting 
halo, whose core radius is too small. We model the Faraday- 
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Figure 17. Profiles of Q/I for the main and counter-jets, inte- 
grated along the jet axis as described in the text. |Q/7| = p for 
an apparent field direction either along or orthogonal to the jet, 
as seen here. Q/I > corresponds to a transverse apparent field, 
Q/I < to a longitudinal one. Full line: main jet; dashed line: 
counter-jet. (a) 69 - 113 arcsec; (b) 113.5 - 157.5 arcsec (the areas 
are the same as in Fig. 171. 



rotating medium as a spherical halo with a core radius 
~225 arcsec and derive an approximate value for the product 
(n /m- 3 ) 2 (B /nT) 2 (Z/kpc) « 700, where no is the central 
density, Bo the central magnetic field and I is the magnetic- 
field correlation length. Our analysis is therefore consistent 
with models of Faraday ro tation proposed for rich clusters 
(e.g. ICarilli fc TavTorll2002l) . but requires much lower densi- 
ties and field strengths. We predict that a tenuous, group- 
scale halo should be detectable in sensitive X-ray observa- 
tions; measurement of its density will allow us to estimate 
the magnetic-field strength. 

We have derived the apparent magnetic field direction 
(corrected for Faraday rotation) and degree of polarization 
at distances between 70 and 160 arcsec from the nucleus. The 
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structure is qualitatively similar to that seen in the flaring 
region, with transverse field on-axis and longitudinal field at 
the edges of both jets, but the degree of polarization on-axis 
is larger. The difference in polarization structure between 
the main and count er-jets observed in the flaring region by 
ICanvin et alJ (120051) persists at larger distances. This can be 
explained fully as an effect of differential aberration on radi- 
ation from intrinsically identical jets, as long as their veloc- 
ities remain significantly relativistic on the relevant scales. 
The asymmetry in RM fluctuation amplitude is consistent 
with the jet orientation required by this analysis and the 
presence of a tenuous, magnetized group halo. 

The large angular size of the flaring region in NGC 315 
and our use of deep observations at several frequencies has 
allowed us to image spectral variations at a level of detail 
not yet achieved in any other jet. Taken together with X- 
ray imaging and modelling of the jet velocity field, this has 
given important insights into the particle acceleration mech- 
anisms. It will be interesting to see whether our results ap- 
ply to other FR I jets and to study the spectral variations in 
NGC 315 over a wider frequency range. 
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